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<Q \ We discuss how to formulate Dirac fermion operator on a finite lattice such 

that it can provide a nonperturbative regularization for massless fermion inter- 
acting with a background gauge field. 
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Consider a massless Dirac fermion interacting with a background gauge field. Our 
present goal is to formulate a nonperturbatively regularized quantum field theory which 
\ at least satisfies the following physical constraints : 



PACS. 11.15.Ha, 11.30.Rd, 11.30.Fs - 
I. Introduction 



(A) In the classical limit, it reproduces the classical physics of the action, 
A = L V>(»7m[^ + igA^(x)}ij(x). 

(B) For topologically trivial gauge backgrounds, and in the weak coupling limit, it agrees 
with the predictions of weak coupling perturbation theory of the action. 



(C) For topologically nontrivial gauge backgrounds, it possesses exact zero modes satis- 
fying the Atiyah-Singer index theorem. 

X. 

Although Wilson's idea H of formulating gauge theories on the spacetime lattice is 
the most successful nonperturbative regularization for pure gauge fields, putting massless 
Dirac fermions on the lattice has been a notorious problem for more than twenty 
years. The resolution of the lattice fermion problem first appeared in the context of the 
Domain- Wall fermion |J, and it motivated the Overlap formalism Q which led to the 
construction of Overlap-Dirac operator || in 1997. We refer to ref. [[J for a recent review 
of the Domain- Wall fermions, and to ref. for a recent review of the Overlap. 

However, if we consider a Weyl fermion interacting with a background gauge field, 
then a completely satisfactory nonperturbative regularization for chiral guage theories ( 
e.g., the standard model ) has not yet been presented up to this moment. 

* Talk given at Chiral'99, "Workshop on Chiral Gauge Theories", Sep. 13-18, 1999, Taipei. 
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In the following, we will concentrate our discussions on the general principles to 
construct chiral Dirac fermion operators on a finite lattice, in particular, for vector-like 
gauge theories such as QCD. With the constraints imposed by the Nielson-Ninomiya no- 
go theorem ||, one can construct a gauge covariant Dirac operator Dona finite lattice 
such that : 

(i) D breaks the chiral symmetry ( i.e., -D75 + 75-D = ) at finite lattice spacing but 

recovers the chiral symmetry in the continuum limit a — > 0. 

(ii) D is local. 

( \D(x, y)\ ~ exp(— \x — y\/l) with I ~ a; or D(x, y) = for \x — y\ > z, where z is 
much less than the size of the lattice. ) 

(iii) In the free fermion limit, D is free of species doublings. 

( The free fermion propagator D^ 1 {p) has only one simple pole at the origin p = 
in the Brillouin zone. ) 

(iv) In the free fermion limit, D has correct continuum behavior. 
( In the limit a — > 0, D{p) ~ i^^p^ around p — 0. ) 

However, one cannot push the property (i) any further, while maintaining properties (ii)- 
(iv). For example, if D is chirally symmetric at finite lattice spacing, then it must violate 
at least one of the three properties (ii)-(iv). We note that these four properties (i)-(iv) 
form the necessary conditions to meet the requirements (A)-(C), however, they are not 
sufficient to guarantee that (C) will be satisfied. 

An example satisfying (i)-(iv) is the standard Wilson-Dirac fermion operator0 

D w = JfJfj, + w (1) 

where 

t»(x, y) = -[U^{x)5 x+ ^ y - Ul{y)5 x ^ v \ , (2) 

o-^t + a u al = 25^ , (3) 



and w is the Wilson term 
a 
2 



w ( x , y) = ~Y. [ 2S ^y ~ U n( x ) S x+(i,y ~ UlWx-fi ,y] ■ ( 5 ) 
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The color, flavor and spinor indices have been suppressed in ([]]). The first term on the 
r.h.s. of (0) is the naive fermion operator which is chirally symmetric at any lattice spacing 
and satisfies properties (ii) and (iv) but violates (iii) since it has 2 d — 1 fermion doubled 
modes. The purpose of the Wilson term w is to give each doubled mode a mass of ~ 1/a 
such that in the continuum limit ( a — > ), each doubled mode becomes infinitely heavy 
and decouples from the fermion propagator. However, the introduction of the Wilson 
term has serious drawbacks. It causes 0(a) artifacts and also leads to the notorious 
problems such as vector current renormalization, additive fermion mass renormalization, 
and mixings between operators in different chiral representations. 

During the last two years, it has become clear that the proper way to break the 
chiral symmetry of D at finite lattice spacing is to conform with the Ginsparg- Wilson 
relation ][| 

D75 + = 2aDR-f 5 D (6) 

where R is a positive definite hermitian operator which is local in the position space and 
trivial in the Dirac space. Then the generalized chiral symmetry (|6]) can ensure that the 



theory is free of above mentioned problems of the Wilson-Dirac operator ||10 



The general solution to the Ginsparg- Wilson relation can be written as [111 



D = D c {l + aRD c )- 1 = (l + aD c R)- 1 D c (7) 

where D c is any chirally symmetric ( -D c 75 + 7s-Dc = ) Dirac operator which must violate 
at least one of the three properties (ii)-(iv) above. Now we must require D c to satisfy 
(iii) and (iv), but violate (ii) ( i.e, D c is nonlocal ), since (|7|) can transform the nonlocal 



D c into a local Dona finite lattice for R = rl with r in the proper range |12|, [TJ|, while 



the properties (iii) and (iv) are preserved. Moreover, the zero modes and the index of 



D c are invariant under the transformation [11, 12]. That is, a zero mode of D c is also a 
zero mode of D and vice versa, hence, 

n+{D c ) = n+{D), n_(D c ) = n_(D), (8) 

index(D c ) = n^(D c ) - n + (D c ) = n-(D) - n + (D) = index(D) . (9) 

Since the massless Dirac fermion operator in continuum is antihermitan, we also 
require that D c is antihermitian ( D\ = —D c ) even at finite lattice spacing. Then 
the chiral symmetry of D c together with its antihermiticity implies that D c satisfies the 
75-hermiticity 

Dl = l5 D cl5 . (10) 
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This implies that D in the general solution (|7]) also satisfies the 75-hermiticity 

D ] = 75^75 • (11) 

Then the eigenvalues of D are either real or come in complex conjugate pairs. Further- 
more, from (0), since R is a positive definite hermitian operator, the lower bound of real 
eigenvalues of D is zero, thus det(D) is real and nonnegative, and is amenable to Hybrid 
Monte Carlo simulation for dynamical fermions with any number of fermion flavors. 

In particular, for R = rl with r > 0, the analysis [ Eqs. (24)-(41) ] in ref. |13|] goes 
through with trivial modifications. The main results are : 

(a) The eigenvalues of D fall on a circle with center at l/2r, and radius l/2r, and have 
the reflection symmetry with respect to the real axis. 

((3) The real eigenmodes ( if any ) at and 1/r have definite chirality +1 or —1. 

(7) The chirality of any complex eigenmodes is zero. 

(5) Total chirality of all eigenmodes must vanish. 
Tr( 75 ) = E s 4>hz<t>s = n+ - n_ + N+ - iV_ = 

where n + (nJ) denotes the number of zero modes of positive ( negative ) chirality, and 
N+(NJ) the number of 1/r modes of positive ( negative ) chirality. From (5), we im- 
mediately see that any zero mode must be accompanied by a real 1/r mode with opposite 
chirality, and the index of D is 

index(D) = n_ - n + = -(JV_ - N+) . (12) 

Now the central problem is to construct the chirally symmetric D c which is nonlocal, 
and satisfies (iii), (iv), and (|10|). Furthermore we also require that D c is topologically 
proper ( i.e., satisfying the Atiyah-Singer index theorem ) for any prescribed smooth 
gauge background. These constitute the necessary requirements jL^] for D c to enter (0) 
such that D could provide a nonperturbative regularization for a massless Dirac fermion 
interacting with a background gauge field. Explicitly, these necessary requirements are : 

(a) D c is antihermitian ( hence 75-hermitian ) and it agrees with 7^(9^ + igA^) in the 

classical continuum limit. 

(b) D c is free of species doubling. 

(c) D c is nonlocal. 
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(d) D c is well defined in topologically trivial background gauge field. 



(e) D c has zero modes as well as simple poles in topologically non-trivial background 
gauge fields ( each zero mode of D c must be accompanied by a simple pole of D c ). 
Furthermore, the zero modes of D c satisfy the Atiyah-Singer index theorem for any 
prescribed smooth gauge background. 

The general solution of D c satisfying these requirements had been investigated in ref. 



18[ and was reported by the author at Chiral'99. However, in general, given any lattice 
Dirac operator D, there exists a transformation T(R C ) for D such that the transformed 
Dirac operator D c = T(R C )[D] is chirally symmetric ||T4| , 15]. Therefore, for the sake of 



completeness, the transformation T(R C ) is outlined in section 2. Then the construction 
in ref. Jl8| is reviewed in section 3. Concluding remarks are briefly outlined in section 4. 



II. A transformation for lattice Dirac operators 

Given any lattice Dirac operator D, in general, there are many different ways to 
construct a chirally symmetric D c out of D. For example, we can construct 

D s = ^{D- l5 D l5 ) (13) 

which is chirally symmetric ( -D s 7s + 75-D s = ). However, this transformation does not 
necessarily preserve the property (iii). For example, if we apply this transformation to 
the Wilson-Dirac operator (|I]), we obtain D s = j^t^, the naive fermion operator which 
suffers from the species doublings. Although D w is free of species doublings in the con- 



tinuum limit, the transformation (|l^) cannot preserve this property since D s is local. 
Therefore, we need a transformaton which preserves the properties (iii), (iv), (|11|) and 
(H), but exchanges the locality of D for its chiral symmetry at finite lattice spacing. The 



transformation [14], [L5| 

T(R) : D ^ D' = T(R) [D] = D(l + RDy 1 = (I + DR)~ 1 D (14) 
which generalizes (0) can serve our purposes. Evidently, the set of these transformations, 



{T(i?)}, form an abelian group with group parameter space {R} |Tj, [15 



For any D satisfying 75-hermiticity, (|TT|), there exists a hermitian R, 



R c = -^(D- 1 + l5 D- 1 l5 ) (15) 
such that 

D c = T(R C )[D] = 2 l5 D( l5 D - D lh )~ l D (16) 
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is chirally symmetric ( 75 -D c + -D c 7s — ) and antihermitian ( D\ = —D c ). 

It is evident that for any two lattice Dirac operators and satisfying (|TT|), 
their corresponding chiral limits obtained from (|T6|), say, and are in general 

different. However, they are related by the transformation 

^TiD^T} (17) 

i 

where each commutes with 75. In general, given any two£>« and D<?\ it is a nontrivial 
task to obtain all Tj in (|17|) . 

For the Wilson-Dirac operator (|IJ), the transformation ( |16D gives 

D c = - w (7^) _1 w • ( 18 ) 

It is antihermitian, nonlocal ( due to the second term ), and satisfies (i), (iii) and (iv). 
Then we can substitute flTBp into (171) with R = rl to obtain a GW Dirac operator 



D 

where 



C(l + r^C)" 1 r CC+CI + ^CC^)- 1 



(19) 



c^^tg-^g- 1 ^. (20) 

On a finite lattice, the GW Dirac operator (|19"D can be constructed to be local but 



not highly peaked in the diagonal elements if the value of r is within a proper range WE 
In other words, if r is too small, then D must be nonlocal since D is close to D c . On the 
other hand, if r is too large, then D becomes highly peaked in the diagonal elements. In 
these two extreme cases, D does not respond properly to the background gauge field ( 
e.g., the anomaly function of D does not agree with the Chern-Pontryagin density of the 
prescribed gauge background ). In the limit the number of sites in each dimension goes 
to infinity, the proper range of r extends to (0, 00). 
The anomaly function of D ( |T9"D can be written 

A D {x) = 2 tr[(l + ^CC^Y 1 - (I + r 2 ^)- 1 ]^, x) (21) 

where tr denotes the trace over the color, flavor and spinor space. 

It is instructive to compare the anomaly function of the Wilson-Dirac operator D w 
(HD to that of the GW Wilson-Dirac operator D (|19D, in a topologically trivial background 
gauge field, on a finite lattice, as shown in Fig. 1 in ref. |~7| and Fig. 1 in ref. |15| , 



respectively. The anomaly function of D w is very different from the Chern-Pontryagin 
density p(x) on a 12 x 12 torus, while that of the GW Dirac- Wilson operator D is in good 
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agreement with p(x) at each site. This demonstrates that the transformation T(r + R c ) 
plays the important role in converting D w into a GW Dirac operator D = T{r + R C )[D W ] 
which is free of 0(a) lattice artifacts, thus D can reproduce the correct anomaly function 
even on a finite lattice provided that the local fluctuations of the background gauge field 
are not too violent. 

However, a lattice Dirac operator satisfying properties (i)-(iv), (|IT]) and the GW 
relation (HD does not guarantee that it has the correct anomaly function in topologically 
nontrivial gauge backgrounds. This can be seen as follows. It is well known that the sum 
of the axial anomaly over all sites of a finite lattice is equal to the index of the lattice 
Dirac operator. Therefore the necessary condition for a lattice Dirac operator to have the 
correct anomaly function is that it possesses exact zero modes satisfying the Atiyah-Singer 
index theorem i.e., n_ — n + = Q. However, not every lattice Dirac operator has exact 
zero modes in topologically nontrivial gauge backgrounds. For example, the Wilson-Dirac 
operator, D w in (|T]), it does not have exact zero modes in topologically nontrivial sectors. 
Then, according to (||), D — T(r + R C )[D W ] also does not have any exact zero modes, 
even though D satisfies the GW relation @, and properties (i)-(iv) and (11). This 



implies that the anomaly function of D w or D = T(r + R c ) [D w ] does not agree with the 
Chern-Pontryagin density for topologically nontrivial gauge fields on any finite lattices. 
Consequently, the disagreement must persist in the continuum limit a — > 0. 

This provides an example to illustrate that any lattice Dirac operator D must possess 
a nonperturbative attribute, the topological characteristics, c[D] |T^, [T7]]. In general, c[D] 
is a rational number, a functional of D and the gauge configuration. In a gauge background 
with nonzero topological charge Q, the index of a lattice Dirac operator is 

n_-n + = c[D] Q . (22) 



In the case of D w and D = T(r + R C )[D W ] (|T9|), n_ = n + = 0, thus the topological 
characteristics c[D w ] = c[D] = 0. For the vacuum sector with nonzero topological charge 
density ( i.e., Q = but p(x) ^ ), c[D] cannot be defined by ([22]), however, it can be 
defined unambiguously |L7| in the context of the anomaly function. 

We note in passing that for the Overlap-Dirac operator D Q = 1 + V ||, the trans- 
formation (T[q) gives D c = 2(1 + V)(l- V)' 1 , as expected. Substituting D c into © with 
R = rl , we obtain the generalized Overlap-Dirac operator 

D = 2(1 + V)[1-V + 2r(l + V)}- 1 . (23) 
III. A construction of D c 

Suppose that we do not have any lattice Dirac operator to begin with. It is still 



possible for us to construct D c [18| according to the necessary requirements (a)-(e) listed 
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in section 1. Since D c is antihermitian, there exists one to one correspondence between 
D c and a unitary operator V, 



1 + V)(l - V)-\ V=(D C - T)(D C + 1) 



-i 



(24) 



where V also satisfies the 75-hermiticity, 75U75 = V'. Then the unitary operator V can 
be expressed in terms of a hermitian operator h, 



V= l5 h 



1 
-1 



hi h 2 
hi h 3 



hi h 2 
-hi -h? 



(25) 



where h\ = hi and h\ = h 3 . Using the unitarity condition V^V = 1, we have h 2 = 1 
1 h\ + h 2 h\ hih 2 + h 2 h 3 \ 



h 2 

Then we obtain 
D = 



h 2 hi + h 3 h 2 h\h 2 + hi 



1 
1 



hi 
h 2 



D R ' 




(I -h 


D L 




-h^il + h) 


-D ] R . The 


general solution to Eq. (^Bj) 


1 






y/1 + tot 




1 


be 16 




VI + 66t 







- x h. 



(26) 
(27) 

(28) 
(29) 



where e ld is an arbitrary phase, and b is any operator. In the following we will restrict 
9 to zero, and also pick the minus sign for h\. Then the general solution for D c can be 
written in the following form 



Dr 



',-1 



1 



-1 



1 + Vl + 66t 



(30) 
(31) 



Due to the presence of the square root in (|30D and (j3~ll) , £) c is nonlocal for nontrivial 
6, thus the requirement (b) is satisfied. In the following, we outline a construction of b 
which satisfies constraints (a)-(d) but not (e). 

In order to have D c satisfy the constraint (a), we first try 



w 



(32) 
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where i M and are defined in 



and (|3|) respectively, and w is a non-singular hermi- 
tian operator which is trivial in the Dirac space and goes to a constant in the classical 
continuum limitQ Then Eqs. ([!(]) and fl3~l~|) become 



where 



D R 
a-t 



(a-ty 
-D\ 



w 



Vi + w- 1 1 2 w- 1 



(33) 
(34) 

(35) 
(36) 



However, t M suffers from species doublings. Now our task is to construct a hermitian 
operator w such that the doubled modes are completely decoupled from the fermion 
propagator D~ l , even at finite lattice spacing. 

As discussed in ref. |18|], we want to construct a hermitian w such that in the free 
fermion limit, it satisfies the following condition, 



w{p) 



> for the primary mode at p = 

< for the doubled modes at p G £$^{0, n/a}\{p = 0}. 



is 



A simplest solution to 

w(p) = c-2^sin 2 (p M a/2), c G (0, 2) 



(37) 



(38) 



Note that the role of w in the general solution of D c is quite different from the Wilson term 
(||) in the Wilson-Dirac operator ([!]). In the general solution of D c , the chiral symmetry 
is always preserved, and the role ofw is to suppress the doubled modes completely at finite 
lattice spacing, while in the Wilson-Dirac operator, the Wilson term breaks the chiral 
symmetry explicitly and gives a mass of order a -1 to the doubled modes such that they 
can be decoupled in the continuum limit ( a — > ). After the gauge links are restored, w 
in the position space becomes 



w(x,y) =c5 : 



x+fi,y 



This is one of the simplest solution of w satisfying the requirement 
limit. Certainly, there exists other solutions to (|37|). 



(39) 



in the free fermion 



t The w operator in this section is different from the Wilson term in Eq. 
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Substituting the D c [ Eqs. (j55|)-(PID with in (Q), and w in with c = 1, and 
the normalization constant 1/2 ] into (0) with R = rl, we obtain the GW Dirac operator 



D 



„2\-l 



B^BB^ + r 2 )- 1 



-B(B^B + r 
r(B^B + r 



a) -i 

2\-l 



where 
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~0 



-1 



The anomaly function of D is 
= 2r 2 tr + 



(StS + r 2 )- 1 



(40) 



(41) 



(42) 



where tr stands for the trace over the color, flavor and spinor space. 

Since the D c of (SO) in the free fermion limit is free of species doubling and has the 
correct continuum behavior, the perturbation calculation in ref. [19| showed that D = 
D c (l + rDc) -1 has the correct chiral anomaly in topologically trivial gauge backgrounds. 
This has been verified explicitly on a finite lattice. An example is shown in Fig. 2 of ref. 

& 



However, for topologically nontrivial sectors, the anomaly function of D also depends 
on the topological characteristics, c[D], which is a non-perturbative attribute of D. Since 
D does not have any exact zero modes in nontrivial sectors, the anomaly function of D 
does not agree with the Chern-Pontryagin density for nontrivial gauge backgrounds. Now 
we have two examples of GW Dirac operators, (B5) and (|4*0"D, which have the correct axial 
anomaly in the vacuum sector but not in the nontrivial sectors. 

In general, if D c is well defined in nontrivial sectors ( i.e., not satisfying (e) ), then 
index [D c ] = index [D] 
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|T2"f ( i.e., c[D] = c[D c ] = 0, D and D c are topologically 
trivial ). On the other hand, if D c has zeros and poles, then c[D] ^ 0, however, it does 
not necessarily imply that c[D] = 1. Sometimes c[D] may even become a fraction ( see 
Table 4 in ref. [[14]] ). The important point is that in general we cannot predict the value 
of c[D] in the topologically nontrivial sectors solely based on the properties of D in the 
free fermion limit or in the vacuum sector. 



IV. Summary and discussions 



The necessary conditions for a lattice Dirac operator D to provide a nonpertur- 
bative regularization for massless fermion interacting with a background gauge field are 
(A)-(C) listed in section 1. The conditions (A) and (B) can be satisfied if D fulfils 
(i)-(iv) and (0), and the chiral symmetry is broken according to the Ginsparg- Wilson 
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relation (||). Since any GW Dirac operator can be represented by the transformation ([/]) 
[ preserving (iii), (iv), (11) and (|j) ] on a chirally symmetric D c , then the conditions for 
D become the necessary requirements (a)-(e) for D c . It should be emphasized that the 
requirement (C) for D ( or the requirement (e) for D c ) cannot be expressed in terms 
of any conditions in the free fermion limit or in the vacuum sector. This naturally leads 
to the concept of topological characteristics [12, 17] associated with each lattice Dirac 



operator D. Although we have no problems to construct D to satisfy (i)-(iv), ([II]) and 
(H), it remains an interesting question whether one can construct a topologically proper 
Dona four dimensional lattice without square root operations. 

If one has a lattice Dirac operator D satisfying (i)-(iv), (|TT| ) and (C), then one can 
transform D into a GW Dirac operator through the transformation D' = T(r + R C )[D] 
which not only preserves all essential physics of D but also can restrict the quantum 
corrections to behave properly such that D is free of 0(a) lattice artifacts, additive mass 
renormalization, etc. 

Recently it has been shown that the effective four-dimensional action ( N s —>■ oo ) 
of the Domain Wall fermion is local and satisfying the GW relation, for gauge fields with 
small field strength [p0[|. However, it is possible to transform the Domain- Wall fermion 



operator into a local GW Dirac operator even at finite N s |LJ| , then the anomalous effects 
( e.g., the residual pion mass in Domain- Wall QCD ) due to chiral symmetry violations 
at finite N s could be suppressed even at moderate N s . Further investigations are required 
before one can tell whether this scenario can be realized or not. 
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